function retval = H_coupling (input_Bo)
 
    global b d p g Rf OA         %%% 常数         
    global Bo II u ii            %%% 时间向量，方位角向量     
    global H W h Ds Dm Os Om     %%% 参数向量          
    global PointSet_three_Bo_u   %%% 三相线点集 
    global PointSet_down_Bo_u_O  %%% 气液面下部点集    
    global PointSet_up_Bo_u_O    %%% 气液面上部点集        
    %%% 特殊 全局变量
    global Point_r_three_Bo_u  Point_z_three_Bo_u  a_down_Bo_u      
    
 %         p*g*H.^3
 %    Bo=------------
 %           b*d
 
 retval=((input_Bo*b*d*10.^6)/(p*g)).^(1/3);
 
end
